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This paper outlines a method for improving the precision of atmospheric neutrino oscillation measurements. 
One experimental signature for these oscillations is an observed deficit in the rate of charged-current 
interactions with an oscillatory dependence on L v / 'E v , where L v is the neutrino propagation distance, and 
E v is the neutrino energy. For contained-vertex atmospheric neutrino interactions, the L v j ' E v resolution 
varies significantly from event to event. The precision of the oscillation measurement can be improved by 
incorporating information on L„/iS„ resolution into the oscillation analysis. In the analysis presented here, 
a Bayesian technique is used to estimate the L v jE v resolution of observed atmospheric neutrinos on an 
event-by-event basis. By separating the events into bins of L v jE v resolution in the oscillation analysis, a 
significant improvement in oscillation sensitivity can be achieved. 
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1. Introduction 

It has now been firmly established by experi- 
ment that atmospheric neutrinos undergo oscilla- 
tions between flavours. In the standard theory, the 
oscillations arise from quantum mechanical mixing 
between the neutrino flavour eigenstates. This is 
governed by a unitary PMNS mixing matrix [U-HI, 
which can be parameterised using three mixing an- 
gles and a CP-violating phase. The amplitudes of 
the oscillations depend on the sizes of the mixing 
angles, whereas the wavelengths depend on the neu- 
trino squared-mass splittings, Am| = m 2 — mf , 
and also on L v jE v , where L v is the neutrino prop- 
agation distance and E v is the neutrino energy. 

The observed atmospheric neutrino data are 
presently well-described by an effective two-flavour 
model of vacuum oscillations. In this approxima- 
tion, the survival probability of an initial state \v^) 
is given by: 



P{Vn ->■ Vy) = l-sin 2 26»f 
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The two-flavour oscillation parameters, Am 2 and 
sin 2 29, have been measured in atmospheric neutri- 
nos by a number of experiments, including Super- 
KamiokandeM, MACRO Q, Soudan 2 and 
MINOS [9h11I]. In addition, the atmospheric neu- 



trino results have been confirmed by the K2K [12J , 
MINOS [lMl and T2K long-baseline experi- 
ments, using accelerator beams of muon neutrinos. 
The MINOS experiment reports confidence limits 



of 



|A<™2 - .•) QO+0.12 



= (2.32+--) x 10-W 



C.L.) and 



sin 2 26» > 0.90 (90% C.L.) [l5|. 



For atmospheric neutrino experiments, one sig- 
nature of neutrino oscillations is a deficit in the ob- 
served rate of charged-current (CC) interactions 
relative to the prediction without oscillations. The 
size of the deficit varies with L v jE v according to 
the formula given in Eq. [TJ The CC interac- 
tions are identified by the presence of an emitted 
muon, which may also be accompanied by addi- 
tional shower activity produced by the final-state 
hadronic system. The neutrino kinematics can be 
.reconstructed by combining measurements of the 
emitted muon and hadronic system, obtained by 
analysing the observed hits in the detector. This 
yields measured values for the neutrino energy, E„ , 
and zenith angle, Q v . The reconstructed zenith an- 
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gle is converted into a propagation distance, L v , 
using a model of atmospheric neutrino production 
height. 

The precision with which an atmospheric neu- 
trino experiment can measure the oscillation pa- 
rameters is dependent on its resolution of the neu- 
trino energy and direction. Typically, the energy 
and direction of the muon can be measured pre- 
cisely. However, the measurement of the hadronic 
system has a worse resolution. Therefore, the res- 
olution of both E v and 6 V improves with E^/ E Vl 
where E^ is the muon energy. The angular resolu- 
tion also improves with neutrino energy since the 
final-state particles are increasingly aligned with 
the incident neutrino direction. This is particularly 
important when the reconstructed muon direction 
is used to approximate the neutrino direction. The 
overall resolution is dependent on the precise con- 
figuration of the detector geometry. For example, 
the resolution is worse if there are fewer hits in the 
detector, or if the final-state particles are not fully 
contained in the detector. 

The resulting resolution of L v jE v varies signifi- 
cantly across the CC data sample. For the mea- 
surement of E v , the resolution is worse for neutrino 
interactions with higher inelasticity, since the muon 
kinematics are better measured than those of the 
hadronic system. For the measurement of L v , the 
resolution is worse at lower neutrino energies, where 
the average scattering angle between the neutrino 
and muon is larger, and also worse at the horizon, 
where L v varies rapidly as a function of zenith an- 
gle, so that a small uncertainty in 9 V translates into 
a significant uncertainty in L v . 

The oscillation sensitivity can be improved by 
incorporating information on the L v jE v resolu- 
tion into the oscillation analysis. An example of 
this technique is the Super-Kamiokande L/E anal- 
ysis [4]. which uses a Monte Carlo simulation to cal- 
culate the average L v jE v resolution on a 2D grid 
of reconstructed neutrino energy and zenith angle. 
A sample of high resolution events is then selected 
by cutting out the regions of energy and angle with 
an average resolution less than 70%. The result- 
ing L„/E v distribution of selected events is seen to 
exhibit a characteristic neutrino oscillation dip. 

This paper describes a Bayesian technique for 
estimating the L u /E u resolution of atmospheric 
and v ^ neutrinos on an event-by-event basis. 
For each event, a probability distribution function 
(PDF) in \og 10 (L u / E„) is calculated by combining 
the measured properties of the emitted muon and 



hadronic system in the event with a Monte Carlo 
simulation of the atmospheric neutrino spectrum, 
neutrino interaction kinematics, and detector reso- 
lution. The L v /E l , resolution is then taken as the 
RMS of this PDF. In the oscillation analysis, all 
selected atmospheric neutrino events are included, 
with events binned according to their resolution. 

The Bayesian technique is demonstrated using a 
simulated atmospheric neutrino data set from the 
MINOS experiment. For the MINOS atmospheric 
neutrino analysis, the separation of events into bins 
of L v j E v resolution is found to yield a significant 
improvement in the oscillation sensitivity. 

2. MINOS Simulation and Reconstruction 



The MINOS Far Detector |17j is a 5.4kton track- 
ing calorimeter located 705 m underground in the 
Soudan mine, MN, USA. Its large mass and under- 
ground location enable MINOS to measure atmo- 
spheric and 77^ disappearance due to neutrino 
oscillations. Atmospheric neutrino and CC 
interactions in the MINOS detector are identified 
by the presence of a muon track with a contained- 
vertex or an upward-going trajectory. The detector 
is also magnetised, allowing muon charge-sign to be 
determined from track curvature. This information 
is used to distinguish between + N — > fjT + X 
and F M + N — » fi + + X CC interactions. 

The analysis presented here uses simulated 
contained-vertex atmospheric neutrino interactions 
in the MINOS detector. The MINOS Monte Carlo 
simulation uses the Bartol 3D |l8[ calculation of 
atmospheric neutrino fluxes and the neugen 
model of neutrino cross-sections. A GEANT3 [2(| de- 
tector simulation tracks the final-state particles and 
provides a full description of the detector response 
and readout. The simulation is used to generate 
a sample of contained-vertex atmospheric neutrino 
interactions corresponding to a total exposure of 
193,000 kton-years. 

A series of dedicated algorithms are used to re- 
construct the muon tracks and hadronic showers 
on the assumption that they are produced by atmo- 
spheric neutrino interactions (2lj . For muon tracks, 
the propagation direction of the muon along the 
track is first determined using timing information. 
A Kalman Filter algorithm is then used to deter- 
mine the muon trajectory through the detector 22j . 
For fully contained muons, which stop in the detec- 
tor, the muon momentum is calculated from the 
track range; for partially contained muons, which 
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exit the detector, the momentum is calculated from 
the track curvature. For reconstructed showers, the 
total hadronic energy is calculated from the visible 
energy in the shower. 

The atmospheric neutrino energy and direction 
are calculated from the reconstructed muon track 
and hadronic shower. The neutrino energy is taken 
to be the sum of the muon and shower energy; and 
the neutrino direction is taken to be the muon di- 
rection. The neutrino propagation distance is cal- 
culated by assuming a fixed 15 km neutrino produc- 
tion height in the atmosphere. 

A set of selection requirements are applied that 
identify contained- vertex muon tracks produced by 
atmospheric neutrino interactions [23j. This yields 
a sample of 5.2 million simulated events which are 
used for this analysis. The and CC compo- 
nent forms 94% of the selected sample of events, 
with the remaining 6% composed of the NC and 
i> e +V e CC backgrounds, which do not oscillate in 
the two-flavour model. The predicted event rates in 
the MINOS detector are 27 events per kton-year in 
the absence of oscillations, and 19 events per kton- 
year using representative oscillation parameters of 
|Am 2 | = 2.32 x l(T 3 eV 2 and sin 2 20 = 1.0, which 
are assumed throughout this paper. 

Fig. [T] shows the true and reconstructed 
\og 10 (L„/ E v ) distributions, plotted with and with- 
out oscillations, for those atmospheric neutrinos 
that pass the selection cuts. For these plots, and 
throughout this paper, the quantities L v and E v are 
measured in units of km and GeV, respectively. In 
the true oscillated distribution, there is a clear dip 
at log 10 (L„ / E v ) w 2.7, corresponding to the initial 
oscillation maximum, and a clear oscillatory struc- 
ture at higher values of log 10 (L J/ / 'E v ), In the recon- 
structed distribution, the oscillations are smeared 
by detector resolution, with only the first oscilla- 
tion dip clearly visible. 

3. Bayesian L/E Analysis 

The aim of the Bayesian analysis is to calculate 
a PDF in \og 10 (L/E) for each event, based on its 
observed kinematics. For MINOS atmospheric neu- 
trinos, these kinematic observables are: the recon- 
structed muon energy, Eff, and muon direction, p^; 
and the reconstructed shower energy, Ef hw . Each 
event is also tagged as q^} = as determined 

by the reconstructed muon charge-sign. The corre- 
sponding true values for these variables are denoted 
Ep, pp, E shw and q„ respectively. 



The required PDF, P, can be written as follows: 

P = P (ft,, E u , q v | p% Ef}, E? hw ,q*) , (2) 

where p v is the neutrino direction, which can be 
converted into the propagation distance, L u . 

Using Bayes' theorem, with the normalisation 
P{p*,E*,E* hw ,q*) = 1, Eq. [Ucan be re- written 
in the following way: 

P = P (p« E* Ef hw: q* | p v , E v ,q v ) 
x P(p v ,E Vl q v ) , 

where the term P(p u , E v , q v ) is the Bayesian prior, 
and gives the expected distributions of p v and E v 
for neutrinos and antineutrinos. 

Eq.[3]is further modified by introducing PDFs re- 
lating the true and reconstructed event kinematics. 
The result is as follows: 

P = P (pf, E*, E? hw ,q* | p M E^E shw ,q v ) 
x P (jpn, E^,E shw | p v ,E v ,q v ) (4) 
x P(p v ,E u ,q u ) . 

In Eq. 31 the first term, which is written 
as P(p$,E*,Ef t hw ,qX'\p li ,E l j l ,E 8 hw,qv), contains a 
set of resolution functions connecting the recon- 
structed muon and shower kinematics to the un- 
derlying true muon and shower distributions. This 
term also contains the relative selection efficiency 
for atmospheric neutrinos as a function of energy. 
The second term, P(p ll ,E li ,E shw \p v ,E y ,q y ), con- 
tains a set of kinematic distributions connecting the 
true muon and shower kinematics to the underlying 
neutrino interaction kinematics. This term incorpo- 
rates the relative cross-sections and kinematic dis- 
tributions for quasi-elastic (QE), resonance (RES), 
and deep-inelastic (DIS) and CC interactions. 
The third term, P{p v , E u ,q v ), is the Bayesian prior, 
described above. Taken together, the three terms 
combine the relative probability of an atmospheric 
neutrino interaction at a given energy and angle, 
with the distribution of kinematic observables pro- 
duced by the neutrino interaction. 

The following subsections describe the calcula- 
tion of each of the terms in Eq. 01 along with the 
approximations made in their calculations. 
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Figure 1: The top panels show the true and reconstructed L V IE V distributions for simulated contained- vertex atmospheric 
neutrinos in the MINOS detector. In each case, the solid line shows the prediction without oscillations and the dotted line 
shows the oscillated prediction, using representative oscillation parameters of Am 2 = 2.32 X 10 — 3 eV 2 and sin 2 29 = 1.0. The 
bottom panels show the ratios of the distributions calculated with and without oscillations. The first oscillation maximum is 
visible in the ratio of reconstructed L v /E v distributions as a dip at log 10 (L v /E v ) ~ 2.7. However, at higher values of L^/E^, 
the oscillations are smeared out. (Note: the up-turn at very high L v /E v , where the lowest energy and hence shortest events 
occur, is due to an increased contamination of downward-going neutrinos mis-identified as upward-going neutrinos.) 
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3.1. Resolution Functions 

The resolution functions translate the underly- 
ing distributions of true muon and shower variables 
into the distributions of reconstructed variables. 
To calculate the resolution functions, it is assumed 
that each observable is measured independently of 
the others. In this approximation, the resolution 
functions decouple into single-variable parameteri- 
sations. For the muon and shower energy, these are 
given by P{E*\E^) and P{Ef hw \E shw ,q v ) respec- 
tively, where separate shower resolution functions 
are calculated for neutrinos and antineutrinos. For 
the muon direction, a perfect resolution is assumed. 
This is a reasonable approximation, since the aver- 
age angular resolution of 1 — 2 degrees for recon- 
structed muon tracks is significantly smaller than 
the average angle between the neutrino and emit- 
ted muon. Therefore, the muon direction is fixed 
by setting P(p^ = S(pj} — p^). For the muon 
charge-sign, the set of probabilities P(q^\q v ) are 
approximated using a similar method. The atmo- 
spheric neutrino selection criteria are used to iden- 
tify muons that have a well-measured curvature. 
For these events, a perfect charge-sign reconstruc- 
tion is assumed, given by P(q^}\q u ) = S(q^ — q v ). 
The vast majority of the contained-vertex events 
fall into this category. For those remaining events, 
it is assumed that P(q^\q v ) = 1/2. 

With these approximations, the overall resolution 
function can be written as follows: 



r {PfJ, ! ) ^Shw 7 % 



Ptj,, E^, E s hw,qv) 

*P{q?\ q») (5) 
xP{E^\E lt ) xP(E^ hw \E shw ,q u ) 

To calculate P{E^\E^) and P{E* w \E shw ,q v ), 
the 5.2 million Monte Carlo events are used to pop- 
ulate 2D distributions in the reconstructed and true 
energy of muon tracks and hadronic showers. The 
distributions are then parameterised to create res- 
olution functions that return a relative probability 
for each combination of reconstructed and true en- 
ergy. For muon tracks, separate parameterisations 
are constructed for stopping and exiting muons, 
which correspond to fully contained and partially 
contained events, respectively. For hadronic show- 
ers, separate parameterisations are constructed for 
Up CC QE, u M CC RES+DIS, V M CC QE and 77 p 
CC RES+DIS neutrino interactions. 

Fig. [5] shows a parameterisation of the average 
fractional resolutions as a function of true energy, 
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Figure 2: Parameterised energy resolution functions for re- 
constructed atmospheric neutrino and 77 M CC interactions 
in the MINOS detector. The fractional energy resolutions are 
plotted as a function of the true energy for: the muon energy 
calculated from range for fully contained muons (solid line) ; 
the muon energy calculated from curvature for partially con- 
tained muons (dashed line); and the hadronic shower energy 
(dotted line). 



for the muon energy, determined from either range 
or curvature, and for the hadronic shower energy. 
The resolutions are averaged over all the selected 
atmospheric neutrino and v ^ CC interactions. 
For fully contained events, the muon energy is de- 
termined from track range with a typical resolu- 
tion of a few percent. However, the muon en- 
ergy from track curvature and the hadronic shower 
energy both have significantly worse resolutions. 
Therefore, the overall neutrino energy resolution is 
better for fully contained than partially contained 
events, and also better for neutrino interactions 
with smaller inelasticities. 

3.2. Kinematic Distributions 

A set of kinematic distributions are used to cal- 
culate PDFs of the true muon and shower observ- 
ables for a given neutrino energy and direction. To 
construct the kinematic distributions, the NEUGEN 
simulation is used to generate a sample of 1 billion 
Up, and CC interactions with a uniform energy 
spectrum. The simulated interactions are separated 
into neutrinos and antineutrinos and are binned ac- 
cording to their interaction type (QE, RES, DIS). 
In each bin, a 3D PDF in neutrino energy, E v , and 
the kinematic variables, W 2 and y, is populated, 
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Figure 3: Distributions of true neutrino energy for Monte Carlo atmospheric (left) and (right) CC interactions that 
satisfy the event selection cuts. In each plot, the open circles show the distributions of selected Monte Carlo events, and the 
solid line shows the parameterisations used as priors in the Bayesian L/E analysis. At low energies, the spectrum rises sharply 
as the selection efficiency increases rapidly from zero; at higher energies, the spectrum falls away, reflecting the underlying 
atmospheric neutrino energy spectrum. 



where W 2 is the invariant mass squared of the final- 
state hadronic system, and y is the inelasticity of 
the interaction. By neglecting the effects of Fermi 
momentum in the interactions, each combination of 
kinematic variables {E v , W 2 , y) can be mapped on 
to a single muon energy, hadronic energy, and rela- 
tive angle between the neutrino and muon. There- 
fore, for a given neutrino energy and direction, 
PDFs of the true muons and shower observables 
(p^jE^^Eshw) can be calculated by marginalising 
the kinematic distributions over W 2 and y, weight- 
ing each bin of interaction type according to its rel- 
ative probability as a function of neutrino energy. 

3.3. Bayesian Prior 

The Bayesian prior gives the expected distribu- 
tions in energy and angle for and v ^ CC inter- 
actions. For the analysis presented here, the distri- 
butions are obtained from the Monte Carlo simula- 
tion. The expected energy distribution, P(E u ,q u ), 
is found by calculating the product of the flux and 
total cross-section for neutrinos and antineutrinos. 
For the angular distribution, it is assumed that the 
incident flux of atmospheric neutrinos is isotropic. 
This is a reasonable approximation, since the se- 
lected events have an average energy of > 1 GeV. 
In this energy region, the atmospheric and 



flux is approximately uniform as a function of the 
neutrino zenith angle [l8j . 

The energy distributions of and F M CC in- 
teractions is combined with the selection efficiency, 
e R (E U: q u ), which has a strong energy dependence. 
Fig. 13] shows the resulting energy distributions, 
P R {E v ,q v ) = P{E v ,q v ) x e R (E u ,q v ), for selected 
neutrinos and antineutrinos. These distributions 
are parameterised for the Bayesian L v / E v analysis. 
In each case, the distributions rise sharply at low 
energies, as the selection efficiency increases rapidly 
from zero; the distributions then fall away again 
at higher energies, reflecting the underlying atmo- 
spheric neutrino energy spectrum. 

The calculated L v jE v resolutions are found to 
be insensitive to the exact choice of Bayesian prior, 
particularly the best L v jE v resolutions, where only 
a small region in neutrino energy and zenith an- 
gle contributes significantly to the Bayesian PDF. 
When a flat energy spectrum is used as a prior, 
similar results are obtained. 

3.4- Calculating the Bayesian PDF 

For each event, a posterior PDF in \og la (L/E) 
is calculated by multiplying the resolution func- 
tions, kinematic distributions and Bayesian prior, 
and marginalising over the kinematic variables E u , 
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Figure 4: Examples of Baycsian PDFs in log lf) (L„ /E v ), log 10 (L^) and \og 10 (E v ) for two reconstructed atmospheric neutrinos. 
The left panels show a partially contained event with reconstructed energies of Ej? = 2 GeV and E^ w = 1 GeV. The muon 
direction is perfectly horizontal and the charge-sign is negative. The L^/E^ resolution is calculated to be 0\ %(L v /E v ) = 0.87. 
The right panels show a fully contained event with reconstructed energies of E^ = 4 GeV and E^ hw = 0.5 GeV, The muon 
direction is directed 45° upwards, and the charge-sign is positive. The L v / E v resolution is calculated to be &iog(L v /E u ) = 0.07. 
In each plot, the vertical lines indicate the reconstructed values, and the histograms give the calculated PDFs. 
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W 2 , y. At each point in this parameter space, the 
true muon momentum and hadronic energy is cal- 
culated, along with the relative angle between the 
muon and incident neutrino. The true muon direc- 
tion is fixed on the reconstructed muon direction, 
and an additional integral is then performed over 
an angle, 0, which rotates the neutrino around this 
direction. At each point in the integral, the true 
values of L u and E„, and therefore log 10 (L„ 
are calculated. A weight is assigned as given by 
the product of the input PDFs, and the value of 
\og 10 {L v /E v ) is entered into the posterior PDF 
with this weight. The values of log 10 (L„) and 
\og w {E v ) are entered into a separate set of PDFs, 
providing information on the relative resolution of 
L v and E v for each event. After the full PDF in 
log 10 (£„/!?„) has been calculated, its RMS value, 
labelled ct\ os (l v /e v )i is taken as the L v jE v resolu- 
tion of the event. 

3. 5. Examples 

Fig. H] shows examples of PDFs in log 10 (X„), 
log 1Q (E v ) and \og 1Q {L v / E v ) for two reconstructed 
atmospheric neutrinos. The first event is a par- 
tially contained muon with reconstructed energies 
of E R = 2 GeV and E* hw = 1 GcV. The recon- 
structed muon track direction, p R , is perfectly hor- 
izontal, and the reconstructed muon charge-sign is 
negative, q R = —1, implying a CC interaction. 
For this event, the resulting PDF in \og w (L v /E v ) 
has a double-peaked structure, arising from the 
rapid variation in the propagation distance, L v , 
with the neutrino zenith angle close to the horizon. 
The resolutions on L v and E v are calculated to be 
^logCL^) = 0.86 and (T\ os (e„) = 0.11, respectively, 
and the overall resolution is <J\ os (l u /e u ) — 0.87. 
This is an example of a low resolution event, where 
the broad \og w (L v ) distribution makes the domi- 
nant contribution to the overall L v jE v resolution. 

The second event is a fully contained muon that 
has reconstructed energies of E R — 4 GeV and 
Ef hw = 0.5 GeV. The reconstructed muon direc- 
tion is directed upwards at an angle of 45° to the 
horizontal, and the reconstructed muon charge-sign 
is positive, q R = +1, implying a CC interaction. 
Since the muon direction is significantly above the 
horizon, the variations in L v as a function of zenith 
angle are small, giving a narrow PDF in log 10 (L J/ ), 
with <7i og (£ iy ) = 0.06. Since the event is fully con- 
tained, the neutrino energy is also well-measured 
and the PDF in log 10 (£' l/ ) returns a small resolu- 
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Figure 5: Distribution of Bayesian L/E resolution, 
(T log(i/B)i f° r the simulated atmospheric neutrino sample. 
The solid histogram shows the predicted distribution in 
the absence of oscillations; the dotted line shows the pre- 
dicted distribution for representative oscillation parameters 
of Am 2 = 2.32 X 10" 3 eV 2 and sin 2 20 = 1.0. 



tion of cri g(B 1/ ) — 0.04. The two distributions com- 
bine to give a narrow PDF in log 10 (L t/ /i?i / ), with 
an overall L V IE V resolution of <J\ og (L v Ie„) = 0.07. 
This an example of a high resolution event. 

4. Separation of Events by L/E Resolution 

Fig. [S] shows the distribution of ciog(L„/.E„) f° r 
the simulated event sample, plotted with and with- 
out oscillations. There is a substantial spread of 
<T iog(L u /E l ,) values across the event sample, corre- 
sponding to ~ 25% of the overall spread of the 
reconstructed neutrino log 10 (L„/E v ) distribution. 
The lowest values of <Ji ok (l u /e u ) are roughly an or- 
der of magnitude smaller than the highest values. 
The shape of the predicted cri og (£„/ e v ) distribution 
is also approximately independent of oscillations. 

To calculate oscillation sensitivities, the events 
are separated into bins of L^/E^ resolution, using 
the two-bin, four-bin and eight-bin schemes given 
in Table [T] Fig. [S] shows the reconstructed distri- 
butions of rog 10 (L„ / E v ) for the four-bin scheme, 
calculated with and without neutrino oscillations. 
Also shown are the ratios of the distributions with 
and without oscillations for each resolution bin. 
The oscillation dip is seen to become increasingly 
pronounced with improving resolution, and is most 
sharply defined in the bin with best resolution. 
Here, the ratio initially falls with log 10 (L v /E u ), 
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Bin number Two resolution bins Four resolution bins Eight resolution bins 



1 0.00 < CTi og (L/£) < 0.50 0.00 < a log(L/E) < 0.25 0.000 < a log(L/E) < 0.125 

2 0.50 < a los(L/E) < 1.50 0.25 < a los(L/E) < 0.50 0.125 < a log{L/E) < 0.250 

3 - 0.50 < a los{L/E) < 0.75 0.250 < a log{L/E) < 0.375 

4 - 0.75 < a los{L/E) < 1.50 0.375 < a log{L/E) < 0.500 

5 0.500 < (T los{L /E) < 0.625 

6 0.625 < a log{L/E) < 0.750 

7 0.750 < a log{L/E) < 0.950 

8 - - 0.950 < a Xog(L/E) < 1.500 



Table 1: Binning schemes used to separate selected events according to the calculated Ly/Ey resolution. 



reaching a first minimum at the point of maximum 
oscillation probability. As log 10 (£„ / ' E v ) increases, 
a second oscillation dip is clearly visible, before the 
ratio tends to an average value of 1 — \ sin 2 29 = 0.5, 
as the frequency of the oscillations becomes rapid. 
Since the highest resolution bin contains a sample of 
events with a better resolved oscillation structure, 
the separation of events into bins of resolution is 
expected to yield a significant improvement in the 
oscillation sensitivity. 

The shape of the \og w (L u / E v ) distribution in 
each L v jE v resolution bin reflects the underlying 
distribution of neutrino energies and angles. The 
events with worse resolution are associated with 
lower energies and angles closer to the horizon, 
whereas the events with higher resolution are as- 
sociated with higher energies and steeper angles. 
Therefore, the central region of the \og 10 {L v / E v ) 
distribution, associated with horizontal neutrinos, 
is increasingly suppressed in higher resolution bins. 
The oscillation structure is most sharply resolved 
in multi-GeV upward-going events, which typically 
have the best L v /E v resolution. 

5. Neutrino Oscillation Sensitivity Study 

The impact of Ly/Ey resolution binning on the 
oscillation measurement is evaluated by performing 
a maximum likelihood analysis on the reconstructed 
log 10 (L^ / 'E v ) distributions. The projected oscilla- 
tion sensitivity is calculated first without resolution 
binning, and then with resolution binning using the 
two-bin, four-bin and eight-bin schemes. 

For this sensitivity study, the Monte Carlo distri- 
butions are scaled to an equivalent exposure of 37.9 
kton- years, matching the analysis described in [llj |. 
The selected events are binned according to their 



reconstructed muon charge-sign, q—(i>,V,X), corre- 
sponding to neutrinos (^), antineutrinos (p) and 
those events with an ambiguous charge-sign (X). 
Each bin of charge-sign is then separated into the 
required number of L v jE v resolution bins. 

The oscillation sensitivities are calculated 
on a 2D grid in the oscillation parameters 
(Am 2 , sin 2 29). At each grid point, the recon- 
structed \og w {L u / E u ) distributions are calculated 
using the selected Monte Carlo events. In addition, 
a set of predicted distributions are calculated 
using representative oscillation parameters of 
Am 2 = 2.32 x 10" 3 eV 2 and sin 2 29 = 1.0. 
These distributions are treated as simulated 
data, and used to evaluate the following negative 
log-likelihood function: 

— In C = /i — n In [i 

-5ZI>ifc ln (/ifc) (6) 

q i,k \ ' 

This log-likelihood function is divided into the fol- 
lowing terms: 

1. Normalisation: The sums — Tilrifj, rep- 
resent the Poisson probability for observing n 
total events, with a total prediction of fi events. 
The sum is taken over the charge-sign bins 
q=(v,V,X). 

2. Shape Term: The terms J2 q J2i,k n ik ln Uik) 
represent the likelihood functions for each of 
the log 10 (L/£') distributions used in the fit. 
The i-sum is taken over the resolution bins; 
the fc-sum is taken over the \og 10 (L/E) bins. 
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Figure 6: The left panels show the reconstructed \og(L/E) distributions for simulated atmospheric neutrinos, plotted with 
and without oscillations, and separated into the following four bins of resolution: 0<<Ti og (£ /£) <0.25; 0.25<a"i og ( i / B ) <0.5; 
0.5<(Ji og ( i/ / B ) <0.75; 0.75<cri O g(^ / '£;)<1.5. In each panel, the solid line indicates the prediction in the absence of oscillations, 
and the dotted line gives the prediction for representative oscillation parameters of Am 2 = 2.32 X 10 -3 eV 2 and sin 2 28 = 1.0. 
The right panels show the ratios of the Monte Carlo predictions with and without oscillations. The oscillation structure becomes 
increasingly clear for bins of higher resolution, and is most sharply defined in the bin with the highest resolution. 
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Within the sum, mu is the observed number of 
events and fik is the relative probability in the 
i th and k th bins. 

Systematic Error Term: Systematic effects are 
incorporated as nuisance parameters, where 
the shift a.j represents the deviation of the j th 
systematic parameter from its nominal value. 
A set of penalty terms, a|/2a^ , are added, 
where the error o~ aj represents the estimated 
uncertainty in the j systematic parameter. 
A total of 10 systematic effects are included in 
the log-likelihood function, to account for sys- 
tematic uncertainties in the atmospheric neu- 
trino flux and cross-section calculations 



The log-likelihood function is minimised with 
respect to each of the systematic parameters. 

By evaluating the log-likelihood function at each 
grid point, a likelihood surface is constructed in 
(Am 2 , sin 2 20) parameter space. The confidence 
levels (C.L.) on the oscillation parameters are then 
calculated assuming Gaussian statistics, where the 
two-parameter 68% C.L. and 90% C.L. are given 
by the locus of points with log-likelihood values of 
— Aln£ = (1.15, 2.30) relative to the central value 
at the input oscillation parameters. 

Fig.[7]shows the resulting 90% C.L. contours, cal- 
culated without any resolution binning, and for the 
case of two, four and eight bins of resolution. The 
use of resolution binning is found to yield significant 
improvements in oscillation sensitivity, particularly 
for the Am 2 parameter. The sensitivity improves 
with each doubling in the number of resolution bins, 
with the improvements becoming smaller each time. 
A set of single-parameter 90% C.L. are calculated 
for the each of the oscillation parameters and are 
found to improve from 1.4 < |Am 2 |/10- 3 eV 2 < 4.9 
and sin 2 29 > 0.79 without using resolution binning, 
to 1.7 < |Am 2 |/10- 3 eV 2 < 3.2 and sin 2 26 > 0.81 
for the case of eight bins of resolution. 
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Figure 7: The projected 90% confidence limits on the oscil- 
lation parameters Am 2 and sin 2 29, calculated without res- 
olution binning, and for the cases of two, four and eight bins 
of resolution. The sensitivities are calculated by scaling the 
contained-vertex atmospheric neutrino Monte Carlo sample 
to a total exposure of 37.9 kton-ycars. The input oscillation 
parameters are Am 2 = 2.32 X 10 -3 eV 2 and sin 2 29 = 1.0, 
as indicated by the star. 



mospheric neutrino data from the MINOS experi- 
ment. The selected events are binned as a function 
of log 10 and also by L v /E v resolution. The 
resolution binning takes full advantage of high reso- 
lution events, which sharply resolve the oscillations. 
It also allows low resolution events to be included 
in the analysis, which contribute to the oscillation 
sensitivity by providing information on the relative 
rate of upward-going and down-going neutrinos. By 
separating events into bins of L v /E v resolution, a 
significant improvement in the oscillation sensitiv- 
ity can be achieved. 



6. Summary 

For atmospheric neutrino and CC inter- 
actions, the L v jE v resolution is crucially impor- 
tant in determining the sensitivity to oscillations, 
but varies significantly from event to event. This 
paper has described a Bayesian technique for es- 
timating the L v /E v resolution on an event-by- 
event basis, which enables an event sample to 
be separated into bins of resolution. The tech- 
nique has been demonstrated using simulated at- 
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